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ABSTRACT 


Convenient  stability  criteria  are  obtained  for  difference 
approximations  to  initial-boundary  value  problems  associated 
with  the  hyperbolic  system  vl.  =  Au  +  Bu  +  f  in  the  quarter 
plane  x  >  0,  t  >  0.  The  approximations  consist  of  arbitrary 
basic  schemes  and  a  wide  class  of  boundary  conditions.  The 
new  criteria  are  given  in  terms  of  the  outflow  part  of  the 
boundary  conditions  and  are  independent  of  the  basic  scheme. 
The  results  easily  imply  that  a  number  of  well  known  boundary 
treatments,  when  used  in  combination  with  arbitrary  stable 
basic  schemes,  always  maintain  stability.  Consequently,  many 
special  cases  studied  in  recent  literature  are  generalized. 
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£•  la&ggas&SB- 

In  this  paper  we  extend  the  results  of  [2]  to  obtain  easily 
checkable  stability  criteria  for  difference  approximations  of 
initial-boundary  value  problems  associated  with  the  linear  hyper¬ 
bolic  differential  system  =  Aux  +  Bu  +  f  in  the  quarter  plane  * 
x  >  0,  t  >  0.  The  difference  approximations,  introduced  in 
Section  1,  consist  of  arbitrary  basic  schemes  —  explicit  or 
implicit,  dissipative  or  unitary,  two-level  or  multi-level  — 
and  boundary  conditions  of  a  rather  general  type. 

The  first  step  in  our  stability  analysis  is  made  in  Section 
2  where  we  prove  that  the  approximation  is  stable  if  and  only  if 
the  scalar  outflow  components  of  its  principal  part  are  stable. 

This  reduces  the  global  stability  question  to  that  of  a  scalar, 
homogeneous,  outflow  problem  which  thereafter  becomes  the  main 
object  of  the  paper. 

Investigating  the  stability  of  the  reduced  problem  our  main 
results  are  restricted  to  the  case  where  the  boundary  conditions 
are  translatory,  i.e.,  determined  at  all  boundary  points  by  the 
same  coefficients.  Such  boundary  conditions  are  commonly  used 
in  practice;  and  in  particular,  when  the  numerical  boundary 
consists  of  a  single  point  the  boundary  conditions  are  transla¬ 
tory  by  definition. 

The  stability  criteria  for  the  translatory  case,  stated 
without  proof  in  Section  3,  are  given  essentially  in  terms  of 
the  boundary  conditions.  Such  scheme-independent  criteria 
elliminate  the  need  to  analyze  the  Intricate  and  often 


0.2 


complicated  interaction  between  the  basic  scheme  and  the  boundary 
conditions;  hence  providing  convenient  alternatives  to  the  well 
known  stability  criterion  of  Gustafsson,  Kreiss  and  Sundstrom  [3], 
which  is  the  basis  for  our  work. 

As  in  [3] ,  we  assume  that  the  basic  scheme  is  stable  for  the 
pure  Cauchy  problem  and  that  the  approximation  is  solvable.  Under 
these  basic  assumptions  —  which  are  obviously  necessary  for 
stability  --  we  obtain,  for  example,  in  Theroems  3*3  and  3-^>  that 
the  reduced  problem  is  stable  if  the  (translatory)  boundary  condi¬ 
tions  are  solvable  and  satisfy  the  von  Neumann  condition  as  well 
as  an  additional  simple  inequality.  If  the  basic  scheme  is  unitary 
it  is  also  required  that  the  boundary  conditions  be  dissipative. 

Having  the  new  stability  criteria,  we  continue  in  Section  3 
to  study  several  examples.  First,  we  reestablish  the  known  fact 
that  if  the  basic  scheme  is  two- level  and  dissipative,  then  out¬ 
flow  boundary  conditions  determined  by  horizontal  extrapolation 
always  maintain  stability.  Surprisingly,  we  show  that  this  result 
is  false  if  the  basic  scheme  is  of  more  than  two  levels.  Next, 
for  arbitrary  multi-level  dissipative  basic  schemes  we  find  that 
if  the  outlfow  boundary  conditions  are  generated,  for  example,  by 
oblique  extrapolation,  by  the  Box-Scheme,  or  by  the  right-sided 
Euler  scheme,  then  overall  stability  is  assured.  Finally, 
for  basic  schemes  (dissipative  or  unitary)  we  show  that  overall 
stability  holds  if  the  outflow  boundary  conditions  are  determined 
by  the  right-sided  explicit  or  implicit  Euler  schemes.  These 
examples  incorporate  many  special  cases  discussed  in  recent 
literature  [1,2,3,4,5,6,9,10]. 
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In  Sections  4  and  5  we  prove  the  results  stated  in  Section  3. 

It  should  he  pointed  out  that  there  is  no  difficulty  in  ex¬ 
tending  our  stability  criteria  to  cases  with  two  boundaries.  In 
fact,  if  the  corresponding  left  and  right  quarter-plane  problems 
are  stable,  then  by  Theorem  5.4  of  [3],  the  original  two-boundary 
problem  is  stable  as  well. 

Thanks  are  due  to  Bjorn  Engquist  and  Stanley  Osher  for  most 


helpful  discussions 


Consider  the  first  order  hyperbolic  system  of  partial  differential 
equations 

(1.1a)  3ha(x,t)/dt  =  Adu(x,t)/dx  +  Bu(x,t)  +  f(x,t),  x  >  0,  t  >  0, 

where  u(x,t)  =  (u^(x,t),...,u^n)(x,t))  '  is  the  vector  of  unknowns 
(prime  denoting  the  transpose),  f  (x,t)  «  (f  ^(x,t),...,f  ^(x,t))'  is 
a  given  vector,  and  A  and  B  are  fixed  n  x  n  matrices  so  that  A  is 
Hermitian  and  nonsingular.  Without  restriction  we  may  assume  that  A  is 
diagonal  of  the  form 

where  A11  and  A131  are  of  orders  Jxi  and  (n-^)x(n-i),  respec¬ 
tively 

The  solution  of  (l.l)  is  uniquely  determined  if  we  prescribe  initial 
values 

(1.1B)  J}(x,0)  -  |(x),  x  >  0, 


and  boundary  conditions 

(1.1c)  tfto.t)--  SuU(0,t)  +  jg(t),  t  >  0, 


where  S  is  a  fixed  lx( n-l)  matrix,  ^(t)  is  a  given  vector,  and 

U.3)  /-(.W . ..<*>)’.  /.(.(•«» . ,«>• 


is  a  partition  of  u  into  Inflow  and  outflow  unknowns,  respectively, 
corresponding  to  the  partition  of  A. 


1.2 


In  order  to  solve  the  initial-boundary  value  problem  (l.l)  by 
difference  approximations  ve  Introduce  a  mesh  size  h  =  Ax  >  0,  k  s  At  >0, 
such  that  X  »  h/k  =  constant .  Using  the  notation  v^(t)  =  v(vh,t)  we 
approximate  (l.la)  by  a  consistent,  two-sided,  general  multi-step  basic 
scheme  of  the  form 


Q-lX^t  +  k)  s  +  v  *  r,r+l,..., 

(JteO 


(1.4a) 


Q  =  £  A.  IT5,  »  =  v 

<r  j__r  jcr  u  ~ 


V  _  — v+1* 


O’  —  —  1,  .  .  .  ,8  , 


where  the  nxn  matrices  Ajff  are  polynomials  in  A  and  kB,  and  the 
n -vectors  by(t)  depend  smoothly  an  f  (x,t)  and  its  derivatives. 

To  solve  (l.4a)  uniquely,  we  provide  initial  values 

(l-4b)  v,(ffk)  »  ^(ffk),  0,...,s,  v  *  0,1,..., 


where  in  addition  we  must  specify,  at  each  time  step  t  =  crk  >  sk,  boundary 
values  ^(t  +  k),  u  =  0,...,r-l.  The  required  boundary  values  will  be 
determined  by  two  sets  of  boundary  conditions,  the  first  of  which  is  ob¬ 


tained  by  taking  the  last  n  -  f  components  of  general  boundary  conditions 
of  the  form 

+  «  -  £  T^v(t-ak)  4  kd(t), 

**■  H  o*0  u  ^ 

2  E^,  n  ■  0,...,r-l,  a  *  -1, . . . , q  , 


(m) 

J(-l) 


where  the  matrixes  are  polynomials  in  A  and  kB,  the  C 

nonsingular,  and  the  n-vectors  d^(t)  are  functions  of  f(x,t),  g(t) 
and  their  derivatives.  If  we  put 


are 


1.3 


,ii(n)  „ll(u) 


M 


3* 


CH  l(u) 

V 


> 

,3El(n) 

'3° 


.  i, 


in  accordance  with  the  partitions  of  A  and  U  in  (1.2),  (1.3),  this 
set  of  conditions  takes  the  form 

(1.4c)  T^l(4V(t  +  k)  +  vj  (t  +  k) 

«  £  [T^l(^(t-crk)  +  lj^jftt-ck)]  +ke  Jt), 

O*=0  K  UK  K 

^©(u)  _  £  ^(n)^,  amlflLf  n s o,...,r - 1. 

For  the  second  set  of  boundary  conditions  we  use  the  analytic  condition 

(l-1*)  vj(t  +  k)  “  S«Jo  (t  +  k>  +  -6(t  +  k> 

together  with  r  -  1  additional  conditions  of  the  form 
(l.4e)  v*(t  +  k)  =  S[Djl(4)vJ(t  +  k)  +  D^^v^t  +  kJJ+kft^t), 


J=1 


H  “  1»*..»^“1» 


where  the  matrices  £«,) 

and  D 


iJr(n) 

3  —  JJ 

ix(n-i),  respectively  —  are  polynomials  in  the  blocks  (A0**) 


of  orders  £  x  £  and 

dfQTv  -1 


and  kB  ,  a,p  *  I,  I,  of  the  matching  partitions 


.-1 


..IIx-l 
(A  )  0 


0  (A11)-1 


kB  ■=  k  ! 

B11 

A 

BnI 

butl/ 

B11  and 

Bnl, 

and  the  i -vectors 

so  that  D*  *  ^  are  homogeneous  1 

v 

£^(t)  are  again  functions  of  f(x,t)»  ,g(t)  and  their  derivatives. 
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1.4 


We  remark  that  (l.4c,d,e)  can  be  solved  uniquely  for  the  required 
boundary  values  v^t  +  k),  ji  =  0,...,r-l,  in  terms  of  neighboring 
values  of  v,  at  least  for  sufficiently  small  k.  Indeed,  since  B 
introduces  an  0(k)  perturbation  of  the  matrix  coefficients  in  (l.4c,e), 
it  suffices  to  prove  this  statement  for  B  =  0 .  But  then,  using  the 
properties  of  and  D*1^,  it  is  not  hard  to  see  that  = 

^  ~  0  and  that  the  are  noosingularj  hence  (l.4c) 

uniquely  determines  the  vectors  v^  (t  +  k),  ^  •  r-l,...,0  (in  that 
order)  and  substituting  in  (1.4d,e)  we  explicitly  cibtain  y^(t  +  k). 

It  *>  0,...,r - 1. 

We  also  remark  that  while  it  is  a  standard  matter  to  construct 
boundary  conditions  of  the  form  (1.4c)  to  any  degree  of  accuracy,  the 
construction  of  (1.4e)  is  less  obvious;  For  example,  using  (1.1)  we 
find  by  induction  on  J  >  1  that 


(1-5) 

and 

(1.6) 


where  the  operators  and  the  vectors  ^j(x,t),  z^(x,t)  are  given 

by 


L 

x 


A_1  “  B)»  ^j(x»t) 

A  ‘S  +  B»  5j(x»t)  - 


=  i:1 

i«=0 

V 

i«=0 


S  (V  A  &=&*(**> 

^  ^Lj ?  Std-i-i 


How,  if  conditions  of  the  form  (l.4e)  are  required  to  p  order  of 
accuracy,  we  take  a  Taylor  expansion  of  u*(t)  and  use  (1.5)  and  (l.lc) 
to  obtain 


1.5 


(1.7) 


n*(t+k)=f)  uT(0,t  +  k)  +  OO^1) 

4  d=o  d*  a*5 

=  s  [(I-t)JH(0»t  +  k)  -^.(O.t  +  k)]1  +  0(hp+1) 

3=0  3*  3 

P  1  / Su^  (0,t  +  k)  +  g(t+k)\ 

=  S  (i  )3  ~  ’  *  1-Zl(«,t.k) 

S-0  J*  L  V  uH(0,t  +  k)  J  *3  J 

+  0(h0+1), 


-tl 


where  [•  ]x  denotes  the  first  l  components  of  the  enclosed  vectors.  We 


see  that  u^Ct  +  k)  depends  on  time  derivatives  of  u^(0,t  +  k)  which,  using 


(1.6),  may  he  replaced  by  space  derivatives  of  (0,t  +  k)  and  u^COjt  +  k). 


Approximating  these  space  derivatives  by  p -order  accurate  linear  combina¬ 
tions  of  u*(t  +  k),...,uI(t  +  k)  and  u?(t  +  k),...,uI(t  +  k),  respectively. 


<0' ”  "/>•" Wp'w  '  ~J  ~  ~ o ■  "  -~p 

we  finally  obtain  (l.4e)  if  u  is  replaced  by  v  and  terms  of  order 
OCh^1)  are  dropped. 


A  concrete  example  of  a  second  order  accurate  boundary  condition  of 
form  (l.ke)  for  the  special  case  B  =  f  =  0,  is  given  in  [2]. 


The  difference  approximation  is  completely  defined  now  by  (1.4)  and 


ve  wish  to  apply  to  it  the  stability  theory  of  Gustafsson,  Kreiss  and 
Sundstrcb  [3].  Trying  to  fit  our  approximation  into  the  form  discussed 
in  [33  we  realize,  however,  that  while  in  the  present  paper  the  vector  b 
of  the  basic  scheme  (l.4a)  is  a  general  combination  of  f  and  its  de¬ 
rivatives,  in  [33  we  have  b  =  f .  Indeed,  the  general  b  admitted  by  us 

here  is  necessary  if  arbitrary  high  order  approximations  to  (l.la)  are 
* 

desired.  Yet,  it  is  not  hard  to  see  that  this  generalization  does  not 
affect  the  results  of  [33.  We  conclude,  therefore,  that  making  the  same 


theory  of  Gustafsson  et  al.  holds  for  our  case  and  we  raise  the  question 
of  stability  in  the  sense  of  Definition  3*3  of  [3]. 

In  Theorem  2.1  below  we  shall  reduce  the  above  stability  question  to 
that  of  a  scalar  outflow  approximation  with  homogeneous  boundary  conditions. 
To  obtain  this  theorem  we  begin  by  recalling  Lemma  10. 3  of  [3]  which  pro¬ 
vides  a  necessary  and  sufficient  determinental  stability  criterion  given 
entirely  in  terms  of  the  principal  part  of  the  approximation,  i.e.,  the 
part  obtained  by  neglecting  B  and  eliminating  all  inhomogeneity  vectors. 
The  mere  existence  of  such  a  criterion  implies  that  for  stability  purposes 
we  may  study  (1.4)  with  b^t)  =  d^t)  =  e^(t)  =  g(t)  =  =  D*^=0; 

hence,  instead  of  (1.4)  we  may  consider  a  basic  scheme  of  the  form 

For  example,  the  Lax-Wendroff  scheme  [73  fcr  (l.la)  is 
irv(t+k)  =  A_iyv.1(t)+Vv(t)+A1Yv+;L(t)  +  Kbv(t),  A0  =  I  +  kB+£k2B2-X2A2, 
A+1  =  £M  +  X2A2-AXk(ABfBA),  b(x,t)  =  [I  +  ^k(B  +  A  ^  +  ^)]f(x,t)  . 


2.2 


(2.0a) 


<j=0 


v  *  r,r+l, . 


»„  -  £  V3  ’ 

j-“r 


v  s  «Xj  • • • ,s ) 


with  initial  values 


(2*lb)  jyv(ck)  —  vy( ok)  ,  a  -  > 


v  =  0,1, • • • , 


and  boundary  conditions 

JHCu)  vB(t  +  k)  „  n  (t  k) 

-1  '  _ a  ~ 


(2.1c) 

x  ,  p.  =  0,. . 

°  j=0  * 

(2.  Id) 

vj(t  +  k)  -  SvJ  (t  +  k). 

(2.1e) 

vj(t  +  k)  «  2  rIlc(u)v?(t  +  k)> 
u  j=i  3 

where  (2.1a)  is  now  consistent  with 

(2.2) 

bS/bt  «=  A  du/dx. 

i1  e  1»  •  •  .  ,T  ~  1» 


and  the  A.  and  C. 

Jo-  jo- 


31  3T(n) 


are  polynomials  in  A  and  in  A 


33  31 


respectively. 

We  thus  obtain. 


UMiA  2.1.  Approximation  (1.4)  is  stable  if  and  oi 

.  i  *1  t  O  i\  -x-vi  — 


asa 


if  its 


2.3 


according  to  the  partition  of  A  in  (1.2) ,  our  next  step  is  to  split  the 
basic  scheme  (2.1a)  and  the  initial  values  (2.1b)  into 


(2.3a) 

Q-1  j£(t  +  k)-  S^V(t-ok)  , 

(7=0 

oil  -  S  A11*1 
"  J.-r  ^  * 

v  c  r, r+1,.... 

(2.3b) 

^(ok)  =  ^(ok),  cr  =  0,...,s, 

v  *=  0)1) •  •  • ) 

and 

- 

(2.4a) 

+  k)  =  S  Q? 1 Z?  (t  -  <rk), 

<JteO 

j=-r  V  E  * 

V  =  r,r+l,ee.. 

(2.4b) 

(ok)  =  ^  (ok),  0=  0,...,s; 

V  «  0  |  1,  Mt) 

thus  viewing  approximation  (2.1)  as  consisting  of  inflow  and  outflow  parts 
given  by  (2.3)  (2.1d,e)  and  (2.4)  (2.1c),  respectively.  Obviously,  (2.1) 

is  stable  if  and  only  if  both  parts  are. 

We  observe  that  the  outflow  part  (2.4).  (2.1c)  is  self-contained  and 
provides,  via  (2.1d,e),  the  boundary  values  v*(t  +  k),  u  =  0,...,r-l. 

We  may  therefore  consider  (2.1d,e)  as  arbitrary  inhomogeneous  boundary 
values  for  the  inflow  part.  So  by  the  argument  involving  Leona  10.3  of 
[3]  preceding  Lemma  2.1,  we  may  replace  (2.1d,e)  —  without  affecting 
stability  —  by  homogeneous  boundary  values 

(2.5)  +  k)  *=  0,  H  “  0,...,r-l. 

This  gives  us  a  new  self-contained  inflow  part,  (2.3)  (2.3),  whose 


stability  together  Kith  that  of  (2.4)  (2.1c)  is  equivalent  to  the  overall 

stability  of  (2.1). 

Since  the  A.  and  are  diagonal,  we  write 

jo-  JCT 

v  -  ^  v  •  • 

and  split  (2.3)  (2.5)  and  (2.4)  (2.1c)  into  n  scalar  components,  each 

of  the  form 


(2.6a) 


ft  ,vy(t  +  k)  =  2  ftffvv(t  -  ok), 
o*0 


v  =  r,r  +  l,.... 


P 

=  S 

J=-r 


V 


IT5, 


(2.6b)  vy(ok)  =  %(<&),  cr  *  0,...,s,  v  »  0,1,..., 

(2.6c)  T^v  (t  +  k)  =  2  T^v  (t-ok),  U  *  0,...,r  -1, 

a  K  o*0  ^ 

where  (2.6a)  is  consistent  with  a  corresponding  component  of  (2.2), 

(2.7)  du/dt  «  a  du/dx,  a  /  Oj 


and  the  boundary  conditions  (2.6c)  are  either  homogeneous,  i.e.. 


(2.8) 


T 


M 

-l 


1,  T 


,00 


0,  u  *  0,. . . ,r  -  1, 


0,. 


>,D, 


or  are  given  by 

»  _  t  ^ 


(2.9) 


J=0 


CJa'  E“  ’  C0(-l)  a 


“1,  •••»<!» 


depending  on  whether  a  <  0  or  a  >  0,  respectively. 

Since  (2.1)  is  stable  if  and  only  if  (2.3)  (2.5)  and  (2.4)  (2.1c)  are 

stable,  and  since  the  latter  are  stable  if  and  only  if  their  scalar 
components  are,  we  obtain  immediately, 
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LEMMA  2.2.  Approximation  (2.1)  Is  stable  if  and  only  If  the  scalar 
components  of  (2.3)  (2.5)  and  (2.4)  (2.1c),  given  by  (2.6)  (2.8)  and 
(2.6)  (2.9).  are  stable. 

In  Section  4  we  shall  prove : 

LEMMA  2.3.  Die  Inflow  approximation  (2.6)  (2.8)  is  unconditionally 
stable. 

This  lemma  —  due  to  Krelss  [4]  in  the  special  case  when  the  basic 
scheme  is  dissipative,  explicit  and  two- level  —  combined  with  the  pre¬ 
vious  two, finally  yields  the  main  result  of  this  section: 

1HE0HEM  2.1.  Approximation  (1.4)  is  stable  if  and  only  if  the  scalar 
outflow  components  of  its  principal  part  are  stable. 

The  above  discussion  implies  that  from  now  on  we  may  reduce  our 
stability  study  to  scalar  approximations  of  form  (2.6)  with  either  (2.8) 
or  (2.9).  We  thus  conclude  this  section  by  stating  the  basic  assumptions 
of  [3]  relating  to  these  approximations  which  will  hereafter  hold  through¬ 
out  the  paper. 

ASSUMPTION  2.1  (Assumption  3.1,  [3]).  Approximation  (2.6)  is  solvable, 
i.e. ,  there  exists  a  constant  K  >  0  such  that  for  each  w  €  jtg(x)  there 
is  a  unique  solution  y  €  Zg(x)  to 

ft-lyv  *  V  v  ■  r»r  + !»•..» 

T^y^  «  w^,  u  ■  0,...,r-l, 

with 


J 
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!ly||  <  k|M|. 

Here,  lg(x)  Is  the  space  of  all  grid  functions  w  «=  {wy  with 


IMf*  bL^ivj2  <  .. 


ASSUMPTION  2.2  (Assumption  5.1  [3]).  lhe  basic  scheme  (2.6a)  is 
stable  for  the  pure  Cauchy  problem,  -  <*>  <  v  <  ®.  That  is ,  putting 


(2.10) 


we  have: 


Vz)  *  v-D  -  ii 2  ‘a.’ 


J  *=-r,... ,p. 


(i)  The  von  Neumann  condition;  i.e. ,  the  solutions  z(|)  of  the 


equation 


(2.11) 


E  &.(z)e  “  *s  o, 


satisfy  |z(Ol<l  for  all  |||<tt. 

(ii)  These  z(|)  which  lie  on  the  unit  circle,  are  simple  roots  of 

(2.n). 

ASSUMPTION  2.3  (Assumption  5.4  [3 J) •  Ifce  basic  scheme  (2.6a)  is 
either  dissipative,  i.e.,  the  roots  of  (2.11)  satisfy 


(2.12) 


MOl  <1,  0  <  |||  <v; 


or  it  is  unitary,  namely 


(2.13) 


Ml)l  «=  1,  III  <  7T. 


Finally,  for  convenience  only  we  make. 


ASSUMPTION  2.4  (Assumption  5*5  (3)). 

a_r(z),  ap(z)  /  0  for  |z|  >  1 


(2.14) 


The  purpose  of  this  section  is  to  provide  easily  checkable  stability 
criteria  for  outflow  approximations  of  form  (2.6)  (2.9).  In  view  of 
Theorem  2.1  this  is  the  key  to  the  overall  stability  question  of  approxi¬ 
mation  (1.4). 

Our  results  —  stated  below  and  proved  in  Section  5  —  are  essentially 
Independent  of  the  basic  scheme  (2.6a)  and  are  given  solely  in  terms  of 
the  boundary  conditions.  These  results,  however,  do  not  apply  to  general 
boundary  conditions  of  form  (2.6c)  (2.9);  Instead  we  are  concerned  in 
this  section  with  the  translatory  case  where  (2.6c)  (2.9)  are  of  the  form 


(3.1) 


T.l^U  ♦  k)  -  2  T  v  (t  -(rk). 


T°  ■  & c*  eJ  ’  c°C-« 


i  o  , 


0, ...  ,r  —  1  • 


As  mentioned  in  the  introduction,  such  boundary  conditions  are 

widely  used  in  practice  since  the  coefficients  c are  independent 

jo 

of  u  and  all  boundary  values  are  conveniently  determined  by  the 
same  procedure.  Especially,  when  the  numerical  boundary  consists 
of  a  single  grid  point  (r  =  1),  the  computation  at  the  boundary 
is  translatory  by  definition. 

We  associate  new  with  the  boundary  conditions  (3.1)  the  boundary 
characteristic  function" 


(3-2) 

where 

(3.3) 


R(e,k) 


2  C.(e)kJ  , 

d-o  3 


cj(z) 


'J(-l) 


<1 

-  £ 
0teO 


s-^c 


J<r» 


J  *  0,...,n. 


This  allows  us  to  state 


THEOREM  3.1.  Let  the  basic  scheme  (2.6a)  be  two-level  and 


dissipative.  Then  the  outflow  approximation  (2. 6a, to)  (3.1)  Is  stable  If 

(3.4)  B(z,ic)  /  0  V  |z|  >  1,  0  <  |k|  <  1. 

EXAMFIE  3.1  (Krelss  [4,  Theorem  6];  see  also  [1]  and  [2,  Example  (4.5)].) 
Let  the  basic  scheme  (2.6a)  toe  two-level  and  dissipative  ,  and  let  the  bound¬ 
ary  conditions  toe  determined  toy  horizontal  extrapolation  of  order  u  -  1,  i .e. , 

(3.5«)  v^t  +  k)  -  2  (*|)(-l)d+1  Vd(t  +  lc)»  l*"  °»...»r-l. 

J*1 

The  boundary  characteristic  function  —  which  for  one-level  boundary  con¬ 
ditions  is  always  z- independent  —  satisfies 

(3-5b)  B(«c)  =  1  -  E  (“)(-l)3+1  J  «  (l-ic)“  /  0,  0  <  |it |  <  1. 

yi  3 

Bence,  (3*4)  holds  and  toy  Theorem  3.1,  (2.6a,to)  (3.5)  is  stable. 

It  should  be  pointed  out  that  Theorem  3.1  is  generally  false  if  the 
basic  scheme  is  of  more  than  two  levels.  Surprisingly,  even  the  well 
known  result  in  Example  3*1  nay  fh.il  to  hold;  namely,  outflow  dissipative 
multi-level  basic  schemes  (s  > 3)  with  boundary  values  determined  by 
extrapolation  of  type  (3* 5a),  are  not  always  stable.  For  example,  con¬ 
sider  the  3-levelj  5-point  basic  scheme 

(3.6)  vy(t  +  k)  =  [I  -  (E  -  I)2(I  -  E‘1)2]vv(t  -  k) 

+  Xa(E-E"1)vv(t),  0  <  e  <  1,  Aa<l  -e,  v  -  2,3,..., 

with  boundary  values  v^(t  +  k),  4  =  0,1,  determined  by  (3*5a).  As 
shown  in  Section  9  of  [  6  ]  the  basic  scheme  is  dissipative,  and  it  is  not 
hard  to  verify  that  the  rest  of  our  basic  assumptions  are  fulfilled  as  well. 
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Yet,  although  condition  (3.4)  of  Theorem  3.1  is  satisfied  as  exhibited 
by  (3.3b),  we  prove  in  Example  4.1  below  that  approximation  (3.6)  (3. 5a) 
is  unstable. 


Despite  the  above  observation  we  can  strengthen  Theorem  3*1  for 
salti-level  dissipative  basic  she  ernes  as  follows. 


THEOREM  3*2.  Let  the  basic  scheme  (2.6a)  be  dissipative.  Then  the 
outflow  approximation  (2.6a,b)  (3*1)  is  stable  if  (3.4)  holds  and  if 

(3-7)  *(*,«  «  1)  +  0  V  |z|  .  1,  z  /  1. 


Evidently,  Example  (3.6)  (3* 5a)  implies  that  the  additional  condi¬ 
tion  (3*7)  is  essential  for  Theorem  3.2. 

Having  stated  Theorems  3-1,  3>2,  we  see  that  when  the  boundary  con¬ 
ditions  (3.1)  are  not  single-level  (as  in  Examples  3.1),  condition  (3.4) 
may  become  a  cumbersome  inequality  in  two  variables,  z  and  k.  Seeking 
a  convenient  alternative  to  these  theorems,  we  extend  the  range  of  p  in 
(3.1)  to  obtain  the  boundary  scheme. 


(3-8) 


T.ivu(t  +  k)  »  S  V^(t  - erk) , 
o*0 


T«=  E  c  E°, 
i<T 


I*  *  0»1»2, ..., 
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and  In  analogy  to  the  definitions  in  Assumptions  2.1,  2.2,  we  introduce: 


DEFUTmCR  3*1.  The  boundary  scheme  (3*8)  is  said  to  be  solvable 
if  there  exists  a  constant  K  >  0  so  that  for  each  y  e  ig(x)  there  is 
a  unique  solution  w  e  ig(x)  to 

(3*9)  ^ 


3A 

with 

M<xM. 

mEKDECncn  3*2.  The  boundary  scheme  is  said  to  fulfill  the  von  Neumann 
condition  if  the  roots  z(g)  of 

(3.10)  B(z,k  -  e1{)  ■  E  c,(z)eiJ*  =  0 

J-0  3 

satisfy  |z(0|  5  1  for  all  |s |  <  tt. 

We  can  now  state, 

THEOREM  3.3  (1st  Main  Theorem).  Let  the  basic  scheme  (2.6a)  be 
dissipative.  If  (3.7)  holds,  and  if  the  boundary  scheme  (3.8)  is  solvable 
and  satisfies  the  von  Neumann  condition,  then  the  outflow  approximation 
(2.6a,b)  (3.1)  is  stable. 

This  result  is  an  extended  analogue  of  the  main  theorem  (Theorem  2.2) 
of  [2]. 

Useful  sufficient  conditions  for  (3.7)  as  well  as  for  the  solvability 
of  the  boundary  scheme  (3*8),  are  given  in  the  next  two  lenmas. 

H9M&  3.1.  Condition  (3*7)  holds  if  any  of  the  following  is.  satisfied; 
(i)  The  boundary  conditions  (3*1)  sure  two-level  (i.e. ,  q  *  0)  and 
accurate  of  order  zero  at  least. 

(±i)  The  boundary  conditions  are  three- level,  accurate  of  order 
zero  at  least,  and  in  addition  R(z  *=  -1,  k  =  1)  f  0. 

IBM  3-2.  (i)  The  boundary  scheme  (3.8)  is  solvable  if 

*  ;)5)  ^  ^ 


0  <  |  k  |  <  1. 


EXAMPLE  3*2.  (Compare  the  special  cases  [3j6*H31  and  [2,  Example  l].) 
Let  the  basic  scheme  (2.6a)  he  dissipative  and  determine  the  boundary  con¬ 
ditions  by  oblique  extrapolation  of  order  w-1: 

(3-U*)  v^t  +  k)  «  £  (p(-l)J+1v^+J[t-(j-l)k],  u  *  0,...,r^l._ 

The  boundary  characteristic  function  associated  with  (3*  Ha)  is  given  by 
(3* lib)  B(z,k)  *  1  -  £  (tf)(-l)‘3+1  z_JkJ  *  (l-z"1*)", 

3«1  J 

so  obviously  (3.7)  holds.  Further,  the  roots  of  (3.11b)  satisfy 

|z(k  -  e11)!  =  le11!  *  1; 

thus,  the  associated  boundary  scheme  —  which  by  Lemma  3*2(ii)  is  solvable  — 
fulfills  the  von  Neumann  condition,  and  by  Theorem  3.3,  (2.6a,b)  (3.11a) 
la  stable. 


EXAMPLE  3-3-  (Compare  the  special  cases  [3,(6. 3c)], [9  ,(3«^)]  and 
[4,  Example  U].)  Let  the  basic  scheme  (2.6a)  be  dissipative  and  let  the 
boundary  conditions  be  generated  by  the  second  order  accurate  Box-Scheme, 


(3.12) 


v^Ct+k)  +  v^+1(t  +  k)  -  Xafv^t  +  k)  -  v^t  +  k)] 


*  V*)  +  Vl^  +  *»lVl(t)  '  “  “  °» *  *  *  »r  “  I* 

3y  Lama  3.l(i),  (3.7)  is  fulfilled;  and  since 

Bellas)]  -  1  +  Re(s)  +  Xa[l  -  Be(s)]  ^0,  |k|  <  1, 

then  by  Loans  3*2(1)  the  boundary  scheme  is  solvable.  The  boundary  character¬ 


istic  function  is 
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The  characteristic  function  for  (3. 13)  is 

fifejic)  *  1  -  z"2  -  Xa(2Kz_1  -  1  -  z“2) 
and  by  equating  to  zero  ue  find  its  roots , 

z(Kce15)  -  e11  -&^P  ,  b(|)  =y/(Xa)2+e’2iStl-(Xa)2]  j 

hence 

*(*  “  X>  *  JlH  /  ei<P»  0<|q,|<7r, 

and  (3.7)  holds.  In  addition,  since  0  <  Xa  <  1  then  |b(|)|  <  1, 

lei  <  «o 

Isis*. 

and  the  boundary  scheme  satisfies  the  von  Neumann  condition.  Finally, 
since  the  boundary  scheme  is  explicit,  Leama  3.2(11)  implies  solvability, 
and  by  Theorem  3.3  stability  follows. 
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We  remark  that  solvability  of  the  boundary  scheme  is  necessary  for 
Theorem  3.3.  To  see  this,  consider  any  dissipative  basic  scheme  with 
zero-order  accurate  boundary  conditions  of  the  form 


(3.1im)  v  (t  +  k)-ev  +1(t  +  k)  =  v  (t)  -  ®v  +1(t)  ,  0>1,  p,  =  0  r..,r-l. 

By  Lemma  3*l(i)»  (3*7)  is  fulfilled.  Also  the  boundary  characteristic 


function  is 


(3.1irt>) 


r(2,k)  =  (l-O(l-eu)  ; 


hence  its  single  root,  z  *  1,  satisfies  the  von  Neumann  condition.  As 
shown  in  Example  k.2,  however,  the  approximation  is  unstable,  which  is 
explained  by  the  failure  of  the  associated  boundary  scheme  to  be  solvable. 
Indeed,  taking  y  =  0  in  (3.9)  we  find  that  the  grid  function 
w  =  {0~^  w0}“_q  with  arbitrary  wQ,  belongs  to  /g(x)  and  satisfies 
(3.9)i  thus  we  have  neither  the  uniqueness  nor  the  boundedness  of  w 
required  by  Definition  3.1. 

Condition  (3*7)  is  also  necessary  for  Theorem  3*3  as  can  be  shown 
by  taking  (3*6)  with  0  <  Aa  <  ^  and  consistent  boundary  conditions  of 

the  form 

(3.15a)  v  (t  +  k)  a  v  (t  -k)  +  2Aa[v  ,(t-k)  -v  (t-k)],  (j,  =  0,l  . 

As  mentioned  before,  the  basic  scheme  is  dissipative,  and  by  Lemma  3*2(ii) 
the  boundary  scheme  is  solvable.  The  boundary  characteristic  function 


(3* 15b) 


R(z,k)  =  1-  z"2[l  +  2Aa(tt  -  1)]  , 
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so  it  is  not  hard  to  verify  that  the  boundary  scheme  satisfies  the 
von  Neumann  condition  (and  in  fact ,  even  dissipative.)  Yet,  as  demon¬ 
strated  by  Example  4*3  below,  (3.6)  (3.15a)  is  unstable.  The  reason 
Theorem  3-3  does  not  apply  in  this  case  is  that  R(z  =  -1  ,  k  =  1)  =  0, 
i.e.,  (3-7)  is  violated. 

So  far  we  have  treated  in  this  section  the  case  where  the  basic 
scheme  is  dissipative.  For  the  general  case  where  the  basic  scheme  might 
also  be  unitary  we  need, 

DEFINITION  3.3.  The  boundary  scheme  (3.8)  is  said  to  be  dissipative 
if  the  roots  of  equation  (3*10)  satisfy  |z(5)|  <  1  for  0  <  |s|  <  r. 

This  enables  us  to  state, 

THEOREM  3.4  (2nd  Main  Theorem).  Let  the  basic  scheme  (2.6a)  be  dissi¬ 
pative  or  unitary,  let  (3«7)  hold,  and  let  the  boundary  scheme  (3-8)  be 
solvable  and  dissipative.  Then  the  outflow  approximation  (2.6a,b)  (3*1) 

is  stable. 

EXAMPLE  3.5.  (Compare  the  special  cases  [3,  (6.3a)],  [8,  (3*2)1 
and  [2,  Example  2].)  Let  the  basic  sehemc  (2.6a)  be  dissipative  or  unitary 
and  let  the  boundary  conditions  be  generated  by  the  right-sided  explicit 
Euler  Scheme, 

(3JL6)  v^t+k)  *  v^t)  +  Xa[vi+1(t)  -  v^t)],  0  <  Xa  <1,  u  =  0,...,r-l. 

The  boundary  characteristic  function  is  now 

R(z,k)  el-  z-1(l  +  Xa( K  -  1) ] , 
and  since  0  <  Xa  <  1,  its  root  satisfies 
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|z(k  =  ei5)|2  =  (Xa)2  +  (l-Xa)2  +  2Xa(l-Xa)cos  | 

<  (Xa)2  +  (l-Xa)2  +  2Xa(l-Xa)  =  1  ,  0  <  | £ |  <n  ; 

hence  the  corresponding  boundary  scheme  is  dissipative .  Moreover,  since 
(3.16)  is  two- level,  first  order  accurate  and  explicit,  Lemmas  3.1(i) 
and  3-2(ii)  imply  that  (3*7)  holds  and  that  the  boundary  scheme  is 
solvable.  The  hypotheses  of  Theorem  3.4  are  fulfilled  therefore,  and 
approximation  (2.6a,b)  (3.16)  is  stable. 

EXAMPIE  3.6.  (Compare  the  special  cases  [8,  (3.3)]  and  [2,  Example  3].) 
Let  the  basic  scheme  (2.6a)  be  dissipative  or  unitary,  and  define  the 
boundary  conditions  by  the  right-sided,  first  order  accurate,  implicit 
Euler  scheme: 

(3-17)  v  (t  +  k)  -Xa[v ;  ,(t  +  k)  -  y ;  (t  +  k)]  =  v  (t),  n=0,...,r-l  . 

The  characteristic  function  associated  with  (3-17)  is  given  by 

E(z,k)  =  1-Xa(«  -  1)  -  z”1  ; 

so  its  root  satisfies 

|z(k  =  ei5)|2  =  [(Xa)2  +  (1  +  Xa)2  -  2Xa(l  +  Xa)cos  5]"1 

<"[(Xa)2  +  (l  +  Xa)2  -  2Xa(l  +  Xa)]“1  =  1,  0  <  |  £  |  <tt  , 

and  the  boundary  scheme  is  dissipative.  Also,  Lemma  3.l(i)  implies  (3*7); 
and  since 

■.[T^M-l+XaU-MO]  /O  ,  |k  I  <1  , 

then  by  Lemma  3«2(i)  the  boundary  scheme  is  solvable.  Thus,  Theorem  3.4 
applies  and  stability  is  assured. 


In  concluding  this  section  we  claim  that  condition  (3*7)  is 
necessary  for  Theorem  3.4  (as  well  as  for  Theorem  3*3*)  For  example, 
consider  the  ncaadissipative,  3-level  Leap-Frog  scheme 

(3.18)  vy(t  +  k)  =  vv(t-k)  +  ^a[vv+1(t)  -  vv_x(t)],  0  <  Aa  <  v=l,2,. 

where  for  the  boundary  condition  take  (3.15a)  with  4=0.  As  mentioned 
earlier,  the  corresponding  boundary  scheme  is  solvable  and  dissipative 
whereas  (3*7)  is  violated.  In  Example  4.4  we  show  that  approximation 
(3-18)  (3*15a)  is  unstable,  thus  proving  our  claim. 

We  conjecture  that  the  solvability  and  dissipativity  of  the 
boundary  scheme  are  essential  for  Theorem  3*4. 


In  this  section  we  use  the  theory  of  Gustafsson  et  al.[3]  to  obtain , 
in  Theorem  4.2  below,  a  preliminary  stability  criterion  for  approximation 
(2.6).  This  criterion  will  be  a  major  tool  in  proving  Lemma  2.3  and 
Theorems  3 . 1  -  3 .4 . 

Following  [3]  we  associate  with  approximation  (2.6)  the  resolvent 
equation 

(4.1a)  ( Q -  r  z-0"1  Q^wy  =  0  ,  v  =  r,r  +  1,..., 

with  boundary  conditions 


(4.2b) 


(*-=0,...,r  - 1  , 


where  z  /  0  is  a  fixed  arbitrary  complex  number.  This  can  be  written 
as 


(4.2)  G(z)w  =  0  , 


where  G(z)  is  a  linear  bounded  operator  on  i2(x)  defined  by  the  left 
hand  sides  of  (4.1).  We  say  that  z  is  an  eigenvalue  of  approximation 
(2.6)  if  (4.2)  has  a  nontrivial  solution  w  =  w(z)  e  ig(x) .  If  z  is 
not  an  eigenvalue  but  there  exists  a  sequence  {w(z)^}j=1  c  ig(x)  with 
l]w(z)^U  =  1  such  that 

G(z)  w(z)^  j  0  , 


we  call  z  a  generalized  eigenvalue  of  approximation  (2.6). 

Having  these  definitions,  we  restate  the  main  result  of  [3]  in  the 


language  of  [5] : 


values  z  with  |z|  >  1. 

Seeking  a  practical  version  of  Theorem  4.1,  we  first  need  the 
following  characterization  of  solvability. 


I£MMA  4.1  (essentially  Osher,  [8].)  Approximation  (2.6)  is 
solvable  if  and  only  if 

(i)  The  difference  equations 

Q-l  \  =  0  ,  v  =  r,r  +  1,... 

(4.2) 

T  1  W|±  =  0  *  n  =  0,...,r  -  1  , 

have  no  nontrivial  solution  w  e  ^2(x). 

(ii)  The  equation 

(4.3a)  Q.jOO  s  etj ( _i)K J  =0 

has  precisely  rQ  solutions  with  0  <  |kj|  <  1  where 
(4.3b)  rQ  =  max{-j  /  °)  • 

(iii)  )  does  not  vanish  on  the  unit  circle,  l.e. , 

<*_!<*)  /  0  »  M  -  1  • 

Proof.  Conditions  (i)  and  (iii)  coincide  with  Osher's  condi¬ 
tions  (d)  and  (g)  in  [8].  Regarding  (ii),  we  note  that  rQ  >  0 
(or  else  the  basic  scheme  is  unnaturally  shifted  to  the  right.) 


4. 


<!.,(«). «'r°  i  vi/'0 

d=-r0  « 

has  a  single  pole  of  order  rQ  at  the  origin.  Since  by  (iii),  Q_1(k ) 
does  not  vanish  on  the  unit  circle,  we  use  the  Argument  Principle  to 
find  that  there  is  no  change  in  arg[Q,_1(it  =  e1*)]  an  £  varies  from 
-T r  to  tt  if  and  only  if  (ii)  holds.  Thus,  conditions  (i)  -  (iii)  are 
equivalent  to  (d)  (e)  (g)  in  [8] »  and  Theorem  I  of  [8]  completes  the 
proof. 

Recalling  the  functions  a^(z)  in  (2. ID),  we  introduce  now  the 
characteristic  equation  of  the  basic  scheme  (2.6a), 

(4.4)  *  P(z,ic)  s  i)  a.(z)it**  =  0  , 

j=-r  J 

whose  r  +  p  roots  k  j(z)  play  a  central  role  in  determining  the 
eigenvalues  of  approximation  (2.6). 

LEMMA  4.2  (compare  Lemmas  5-1  and  5*2,  [3l  •)  For  |z|  >  1, 
the  characteristic  equation  (4.4)  has  precisely  r  roots  with 
0  <  |kj(z)|  <  1,  p  roots  with  |kj(z)[  >  1,  and  no  roots  with 
|*j(jO|  =  1. 

Proof.  By  Assumption  2.4,  the  leading  coefficients  of  P(z,k) 
do  not  vanish  for  |z|  >1;  hence  (4.4)  has  r  +  p  roots,  all 
satisfying  |kj(z)|  >  0.  Since  these  roots  are  the  solutions  of  the 
the  polynomial  equation 

(4.5)  Pr(z>K)  s  «rP(z,iO  *  0  , 


we  may  study  (4*5)  rather  than  (4.4). 


By  Assumption  2.2  (i). 


Pr(z,K  *  e15)  /  0  ,  |zj  >  1  ,  ||(  <tt  , 

i.e.,  for  |z|  >  1,  P  ( z  ,k  )  does  not  vanish  on  the  unit  circle  |k|  =  1 

Since  the  roots  of  Pr(z,ic)  are  continuous  functions  of  z,  it  follows 

that  for  |z|  >1  the  number  of  roots  satisfying  0  <  |k  (z)|  <1.  is.. 

o 

independent  of  z.  In  particular,  consider  the  limit  case 


’^c)  =  ^  a.,  ,o«J+r 


d=-r  H-1* 

By  Lemma  4.1  (iii), 

Pr(z -»»,  k  =  e15)  ^  0  ,  U|<tt  , 

so  by  continuity  again,  the  number  of  roots  of  Pr(z,n)  satisfying 
0  <  |kj(z)|  <  1  may  be  determining  by  counting  the  roots 
■tj,  |xjJ  <  1,  of  Pr(z  -»»,*)•  We  have 

r-rn  p  3+rn 

Pru  «  s  v-u" 

J  0 

where  rQ  is  defined  in  (4.3b).  Moreover,  by  Lemma  4.1(ii), 
j  J+rQ 

-  r  aj(  -1)K  °  ^nce  £j__r  aj(_i)*  have  precisely  rQ  roots 

with  |k  |  <1.  Thus,  with  its  additional  r  -  rQ  zeroes,  Pr(z-*»,n) 
has  r  roots  with  Ik  I  <  1  and  the  lemma  follows. 


According  to  the  above  lemma  the  roots  of  the  characteristic 
equation  (4.4)  split  for  |z|  >  1  into  two  groups :  r  inner  roots 
satisfying  0  <  |*j(z)|  <  1  and  p  outer  roots  with  |itj(z)|  >  1. 
By  continuity,  therefore,  these  groups  of  inner  and  outer  roots 
remain  well  defined  for  |z|  >1  as  well,  where  milder  inequalities, 


lKj(z)l  and  |*^(z)|>  1,  hold,  respectively*  Since  by  Assumption 
2.4,  k  m  0  is  never  a  root  of  P(z,k)  for  |z|  >1,  we  finally 
obtain: 

I£MMA  4.3.  For  |z|  >  1,  the  r  +  p  roots  tc ^ ( z)  of  the  charac¬ 
teristic  equation  (4.4)  split  into  r  inner  roots  with  0  <  |ti  j(z)|  <  1 
and  p  outer  roots  with  I K  j( z) I  — 

Mow,  let  z  be  given.  It  is  well  known  (e.g.  [5])  that  z  is 
an  eigenvalue  or  a  generlized  eigenvalue  of  approximation  (2.6)  if  and 
only  if  equations  (4.1)  have  a  nontrivial  solution  of  the  form 

C..6)  »v  -  £  Jj,  Taf>  ’  v  -  °-1’2 . 

where  *a(z),  1  <  a  <  N,  are  the  distinct  inner  roots  of  the  character¬ 
istic  equation  (4.4)  each  with  multiplicity  Mq  =  Ma(z).  Here,  $a^(v) 
sure  arbitrary  polynomials  in  v  with  deg[fa^(v)]  =  p,  and  are 
coefficients  whose  number,  by  Lemma  4.3,  is  precisely 


£  M  =  r  . 

0=1 

To  find  the  t  we  substitute  (4.6)  in  (4.1b)  and  obtain  a  linear 

up 

homogeneous  system  of  r  equations  with  r  unknowns. 


(*•7) 


j(z)T'  =  0  , 


where  J(z)  is  an  r  x  r  coefficient  matrix  and  T  =  (^p)  is  the 
unknown  vector.  Obviously,  w  =  w(z)  in  (4.6)  does  not  vanish  if 
and  only  if  (4.7)  bas  a  nontrivial  solution  t;  hence,  z  is  an 
eigenvalue  or  a  generalized  eigenvalue  of  approximation  (3.6)  if  and 
only  if  det  J(z)  /  0. 


This  observation  combined  with  Theorem  3.1  gives  the  following 


equivalent  of  Lemma  10. 3  of  [3]: 


U3MA  4.4  (Gustafsson  et  al.  [3].) 
stable  if  and  only  if 


(2.6)  is 


det  j(z)  /  0  ,  |z|  >  1  . 

How  for  z,  | z |  >  1,  with  corresponding  distinct  inner  roots 
*a  *  Ka(z)»  1  <  a  <  H,  each  with  multiplicity  M^,  we  a 

specific  choice  for  the  polynomials  $a^v)  1“  (4.6) : 


Thus,  (4.6)  becomes 


H  V1 


“v  -  £  £  ^  (S)-« 


,  v  «  0 ,1,2, . . . , 


and  substituting  in  (4.1b)  —  with  T^  given  by  either  (2.8)  or 
(2.9)  -  we  obtain  explicit  expressions  for  the  system  J(z)t'  » 0,  namely 


(4.8) 


N  Maf 1 


a*l  p=0 


E>  P!(p)  “«**<»»  ’  ,‘*0”"’r‘1  • 


(4-9)  l  *£  l  (•&>  -£ *'°'1  to  - » 

according  to  the  inflow  and  outflow  cases  (a<0  or  a>0), 
respectively.  To  further  simplify  these  expressions  we  use  the 
boundary  coefficients  to  introduce  r  boundary  functions 


yz,*)  m 


(4.10) 


a  <  0 
a  >  0  , 


*  *  0,...,r-  1  . 


9 


a  >  C 


a=0 


-0-1 


systems  (4.8)  and  (4.9)  both  take  the  unified  compact  form 


N  a 


r  J 


£  r 

Qt=l  (3=0 


dKR. 


**  1 


^a(z) 


*  hfi 


=  0 


|x  =  0 , .  • .  ,r  *  1 


so  the  coefficient  matrix  J(z)  can  be  conveniently  written  as 
(4.11a)  J(z)  —  [H(z jtt ^>M^)  > .  •  •  jH(z,Kjj,M^) ]  , 


where  H(z,K^,Ma) ,  1  <  a  <  N,  are  the  r  x  blocks 
(4.11b)  H(z,Ka,Ma) 


R0(z,k) 

Bo(z,k) 

R^z,*) 

R^ZjK) 

d 

R^(  z,k) 

M_  - 
d  0,-1 

Rx(z,k) 

• 

• 

•  i 

*  d« 

• 

• 

•  • 

» •  *  *  *  M  , 
d*  0-1 

• 

• 

• 

Rr-1^Z,K) 

Rr-l(z,,0 
•  M 

Rj,.^(  z,k) 

L 


K*Ka 


This  expression  for  J(z)  gives  us  a  concrete  analogue  of 
Lemma  4.4: 

LEMMA  4.5.  Approximation  (2.6)  is  stable  if  and  only  if  f 
every  z,  |z|  >  1, 

(4.12)  det  J(z)  s  det[H(z)K1,M1),...,H(z,KR,MR)]  j*  0 

where  *a  =  *a(z)>  1  <  a  ^  N,  ftgg-frhft-fllatilBrt  iSBSL  ISSSa.  9f 
each  with  multiplicity  ■  *^(*)» 


4.8 


A  milder  version  of  this  theorem  is  given  in  Theorem  3*2,  [2]. 

We  return  now  to  the  case  where  the  inflow  boundary  conditions 
are,  as  before,  homogeneous  and  given  by  (2.6c)  (2.8),  whereas  the 
outflow  conditions  (2.6c)  (2.9)  are  translatory  as  described  in 
(3.1).  In  this  case  —  where  both  the  inflow  and  outflow  conditions 
are  of  translatory  type  —  the  boundary  functions  in  (4.1£>)  become 


a  <  0 
a  >  0  , 
.  ,r  -  1  . 


Hence,  denoting 

R(2,k)  s  Rq(z,k ) 

(note:  for  a  >  0  this  definition  coincides  with  the  one  in  (3-2)), 
we  find  that 

R^**)  -  ^(z,*)  ,  u  =  0,...,r -  1  ; 
so  the  blocks  of  j(z)  in  (4. lib)  became 


(4.14)  H(z,«o,Ma) 


R(*»<0 

kR(z,k) 


■r(z,k) 

‘  R(z,  k) 

kR(z,k) 

0 

M -1 

a  a 

kR(z,k) 

• 

• 

icr-1R(z,it) 

V1 

dK 

• 

• 

*r"1R(z,K) 

»  - 

This  allows  ujb  to  simplify  Lemma  4.5  as  follows. 


THEOREM  4.2  (com 


4.9 


pare  Theorem  4.1,  [2].)  Approximation  (2.6)  with 
translatory  boundary  conditions  given  by  either  (3.1)  for  a  >  0  or 
(2.8)  for  a  <  0,  is  stable  if  and  only  if  for  every  z,  |z(  >  1, 
with  corresponding  inner  roots  *a(z),  1  <  a  <  H,  we  have 

(4.15)  ^(z.Kq)  t  0  ,  a  =  1, . . . ,N  . 

Proof.  Suppose  R(z,*a)  =  0  for  some  z,  |z|  >1,  with  a 
corresponding  inner  root  *a.  Then  clearly,  the  first  column  of 
HU,Ka,Ma)  in  (4.14)  vanishes;  thus  det  J(z)  =  0,  and  by  Lemma  4.5 
we  have  instability. 

Conversely,  let  (4.15)  bold  and  take  an  arbitrary  z,  |z|  >1, 
with  distinct  inner  roots  t*a(z),  1  <  a  <  N.  To  prove  stability  it 
suffices,  by  Lemna  4-5,  to  verify  that  the  rows  of 

(4.16)  J(z)  **  [H(z,k 2_jM^) , •  •  •  ,H(z,Kjj,Mjj)] 


are  linearly  independent,  where  the  H(z,!*a,Ma) 
For  that  purpose,  let 


(4.17) 


are  given  by  (4.14). 


=  0 


be  a  vanishing  linear  combination  of  the  rows  of  (4.16),  and  let  us 
rewrite  (4.17)  as  r  scalar  equations 


8ince  by  hypothesis 


4.10 


R<Z»K>lK-Ka^° 


. 

! 

1  <  a  <  H  , 

"  ~  j 


we  expend  the  partial  derivatives  in  (4.18)  by  Leibnitz'  rule  and 
use  induction  on  j  >  0  to  find  that  the  sum  in  (4.18)  must  have 
vanishing  derivatives)  i.e., 


r-1 

r 

)JL*=0 


=  K . 


a 


0 


l<a<N,  0  <  j  <  M  ^ 


Consequently,  the  polynomial 


r-1 

*(k)  h  S 

JJL  =0 


which  is  of  degree  r-1  at  most,  has  r  roots  (*a,  1  <  a  <  N, 
each  with  multiplicity  M^)  so  ¥(k)  s  o  and  the  coefficients 
7^  must  vanish.  By  (4.17),  therefore,  the  rows  of  (4.16)  are 
linearly  Independent  and  stability  follows. 


The  proof  of  Lemma  2.3  and  the  counterexamples  of  Section  3 
are  almost  at  hand  now. 


Proof  of  Lemma  2.3.  By  (4. ID),  the  boundary  function 
associated  with  the  homogeneous  boundary  conditions  (2.6c)  (2.8)  is 

"  R(*»*)  =  ^(*»K)  “1  • 

Thus,  (4.15)  bolds  trivially,  and  by  the  last  theorem  approximation 
(2.6)  (2.8)  is  stable. 


SXAMP1£  4.1.  Consider  the  dissipative  basic  scheme  (3.6) 
with  the  boundary  conditions  in  (3.5a).  The  boundary  function  is 
given  by  (3.5b)  and  for  i  -  -1  it  can  be  shown  (as  in  Lama  6.2, 


[3])  that  the  characteristic  equation  has  exactly  one  inner  root 
satisfying  k(z  >  -1)  *  1.  Hence,  R(z  =  -1,  it  =  1)  «  0,  and  by 
Theorem  4.2  we  have  instability. 

EXAMPLE  4.2.  Take  the  zero-order  accurate  boundary  conditions 
(3.14a)  in  combination  with  any  basic  scheme  (dissipative  or  even 
unitary.)  By  (3.14b),  R(z  =  1,k)  =  0  for  all  it;  so  at  z  =  1 
the  characteristic  boundary  function  vanishes  for  all  inner  roots, 
and  Theorem  4.2  implies  instability. 

EXAMPLE  4.3-  Take  the  same  basic  scheme  as  in  Example  4.1 
with  the  boundary  conditions  in  (3 -15a) .  As  in  Example  4.1  we  have 
an  inner  root  it(z  *  -1)  =  1  for  which,  by  (3.15b),  R(z,k)  =  0. 
Hence,  (4.15)  is  violated  and  approximation  (3-6)  (3*15)  is 
unstable. 

EXAMPLE  4.4.  Consider  the  Leap-Frog  scheme  (3*18)  with  a 
boundary  condition  as  in  Example  4.3.  In  Lemma  6.2  [3]  it  is 
shown  that  the  characteristic  equation  of  (3*18)  has  a  single  inner 
root  it(z  «  -1)  =1.  So  as  in  the  previous  example,  R(z  =  -1,  k  *=1)  =0, 
and  by  Theorem  4.2  instability  follows. 


5. 


!•  Proof  of  Main  Results . 

We  turn  now  to  prove  the  results  stated  in  Section  3>  beginning 
with  the  following  lemma. 


U34MA.  5*1-  For  z  =  1,  the  characteristic  equation  (4.4)  has 
exactly  one  root  satisfying  it(z  =  1)  =  1.  In  the  outflow  case 
(a  >  0)  this  is  always  an  outer  root. 


Proof.  Since  the  basic  scheme  (2.6a)  is  consistent  with 
au/at  =  adu/ax  ,  a  ^  0  , 

the  coefficients  a  must  satisfy  the  ordinary  consistency  conditions 

So 

$  a  .  .  .  S  $  a 

j=-r  c=0  j=-r  •** 


and 


^  Jad(- 
J=-r 


-x  «  S  ^  da.,  -Aa  S  (o  +  1)  a.  , 
'  c=0  j=-r  00  <7=0  j=-r  00 


which  can  be  written  as 


(5-la) 

£  ajU)l 

(5.1b) 

^  3a.(z) 

d=-r  J 

z»l 


I  "  "Xa  ^  £  *\}(z) 

I  z=l  j*-r 


z«l 


or  equivalently  as 
(5.2a)  P(2,k)| 


z=k  =1 


0  , 


at  p(z»*)|  *  -Xa  ^  p(z>K)j 

I  z=*  «1  I 


(5.2b) 


zwcal 


zwc  *1 


5.2 


Here,  a^(z)  P(z,it)  are  defined  in  (2.1)  and  (4.10,  respectively. 

£y  (5.1a),  z  =  1  is  a  solution  of 

^  a.(z)  =  0 
d=-r  3 

and  by  Assumption  2.2(ii)  this  solution  is  simple.  Hence, 


te  $  aj(z)  , 

j=-r  °  z=l 


t  0 


so  by  (5 .2b) 


(5-3)  P(z,k)|  =  -Xa  ^  P(z,k)| 

I  Z=tt  =1  °  Z=K  =1 


^  A  ai(z) 


Z=1 


Having  (5 .2a)  and  (5*3),  we  employ  the  implicit  function  theorem  to 
find  that  in  the  neighborhood  of  z  =  1  the  characteristic  equation 
(4.4)  can  be  uniquely  solved  for  it  as  a  differentiable  function  of 
z  such  that 

(5*4)  k(z  s  l)  s  1  . 


This  is  the  first  part  of  the  lemma. 

To  complete  the  proof  consider  the  outflow  case  a  >  0.  Then 
(5* 2b)  yields 


(5-5) 


d)t(z) 

dz 


z=l 


_  &p  /dP 

**  dz  /  3it 


Z=K  =1 


=  ~  >  0 
Xa 


so  by  (5*4)  (5*5),  for  z  =  1  +e  with  sufficiently  small  e  >  0, 


it(z)  *  1  +  (Xa)"1  e  +  0(e2)  >  1  . 

That  is,  for  z  in  the  right  real  neighborhood  of  z  *  1, 


|*(z)|  >  1  , 


5-3 


and  by  Lemma  4.2  this  inequality  is  valid  for  all  z,  |z|  >  1.  By 
definition,  therefore,  *(z)  of  (5*4)  is  an  outer  root  of  (4.4)  and  the 
lemma  follows. 

Proof  of  Theorem  3.1.  Take  an  arbitrary  z,  |zj  >1,  and  let 
=  kq(z)  be  s  corresponding  inner  root.  In  order  to  prove 
stability  it  suffices,  by  Theorem  4.2,  to  show  that 

(5.4)  /  0  . 

By  Lemma  4.3»  we  have  0  <  |k  (z)|  <  1,  where  for  0<|*  (z)|  <1, 
(5.4)  is  implied  by  (3-4).  Hence,  we  may  restrict  attention  to  inner 
roots  on  the  unit  circle,  i.e.,  k^(z)  “  e*^,  |||  <  tt.  Since  the  basic 

scheme  is  dissipative  then  by  (2.12),  the  solutions  z  =  z(k)  of  (4.4) 
satisfy 

(5-5)  \z(k  =  ei5)|  <  1  ,  0  <  |e|  <tt  . 

Thus,  for  |z|  >1,  (4.4)  has  no  roots  k  =  e*^,  0  <  |||  <  ir,  and 
our  discussion  is  further  reduced  to  | z)  >  1,  Ka(z)  =  1.  Next,  by 
continuity,  (5-5)  yields 

|z(k  =  1)|  <  1  J 

so  k  «  l  is  ruled  out  as  an  inner  root  for  jz|  >  1  and  it  remains 
to  consider  |z|  =1,  *a(z)  =  1.  Finally,  by  Lemma  5.1,  k  =  1  is 
is  excluded  as  an  inner  root  for  z  =  1  and  we  are  left  with 

(5*6)  |s|  -  1  ,  s  ft  1  ,  kq(z)  =  1  . 

Since  the  basic  scheme  is  consistent  then  by  (5 .2a) , 

(5-7)  P(z  -  1,  k  -  1)  =  0  . 
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Moreover,  since  the  basic  scheme  is  two-level,  P(z,ic)  is  a  poly¬ 
nomial  of  first  degree  in  z"^  where,  by  (5*7)»  its  only  root  is 
z-1  -  1.  Thus, 

P(z,Kal)j/0,  |  z  I  =  1  ,  Z  ^  1  , 

and  the  proof  is  complete. 


For  Theorem  3*2  we  repeat  the  previous  proof  to  the  point  where 
the  remaining  values  of  z  and  *a(z)  to  be  studied  are  given  in 
(5.6).  At  this  point,  (3*7)  implies  (5.4)  and  stability  is  assured. 
To  prove  Theorems  3*3  and  3*4  we  need  yet  another  result. 


IEMMA  5-2.  Let  the  boundary  scheme  (3-8)  be  solvable  and 

satisfy  the  von  Neumann  condition.  Then, 

(5.8)  R(z,k)  /  0  ,  |k|  <  1  ,  |z|  >  1  , 

and 

VT9)  R(z,k)  ^  0  ,  Jk|  <  1  ,  |z|  >  1  . 

Proof.  Apply  Lemma  4.1  to  the  solvable  boundary  scheme,  rather 
than  to  the  basic  scheme,  with  equation  (4.3a)  replaced  by  its 
boundary  counterpart 


m 


T-i(,)  *  £  cj(-dkJ  - 0 


Since  the  boundary  scheme  is  one  sided  where  by  (3-1)  cj(-i)  ^ 
then  parts  (ii)  -  (iii)*  of  the  lemma  imply  that 


We  also  have, 

T.iOO  >  o<  M  <1 

so  all  told 

T-l<°>  **  c0( -1)  *  0  5 

(5.10) 

,  M«£l  . 

A 
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Since  the  boundary  scheme  satisfies  vcn  Neumann  condition  then 
by  (3* ID),  for  Jz|  >1,  R(z,k)  does  not  vanish  on  the  unit  circle 
| k |  =1.  Hence  —  as  in  the  proof  of  Lemma  k.2  —  the  number  of 
roots  of  R(z,k)  satisfying  |n(z)|  <1  is  independent  of  z,  |z|>l, 
and,  by  continuity,  equals  that  of  the  roots  k,  |k|  <  1,  of 

R(z  — » co,  k)  =  T_^(k)  =  0  . 

By  (5.10),  therefore,  R(z,k)  has  no  roots  | »c ( |  <  1  for  ]z|  >  1, 
and  (5*8)  bolds. 

To  obtain  (5.9),  we  merely  note  that  by  (5*8),  the  roots  k(z) 
of  R(z,k)  satisfy  |tt(z)|  >1  if  |z|  >  1.  Thus,  for  |z|  >1, 
those  continuous  roots  satisfy  (tt(z) |  >  1  and  the  lemma  follows. 

Proof  of  Theorem  3-3.  Since  the  boundary  scheme  is  solvable  and 
satisfies  the  von  Neumann  condition.  Lemma  5.1  implies  (3.1+),  and  by 
Theorem  3*2,  approximation  (2.6a,b)  (3>1)  is  stable. 

1 

Proof  of  Theorem  3 .^ •  As  in  the  proof  of  Theorem  3.1,  let  z 
satisfy  |z|  >1,  let  Ka(z)>  lKa(t)|  <  1  ,  be  a  corresponding 
inner  root,  and  let  us  prove  that 

(5-11)  R(z,»ta)  ^  0  . 

By  Theorem  k.2 ,  this  will  Implies  stability. 

Comparing  Definitions  3*2  and  3*3>  we  immediately  see  that  since 
the  boundary  scheme  is  dissipative,  it  satisfies  the  von  Neumann 
condition;  so  Lemma  3*2  applies,  and  it  remains  to  verify  (5*U)  f°r 
z  and  kq(z)  with 

1*1  -  1  >  ha(*)l  “  1  • 


duutfUfiMkMiitlii 


Indeed,  for 


M  =  1  >  lKa(z)i  =  1  »  Ka(z)  i  1  » 

(5-U)  follows  from  the  dissipativity  of  the  boundary  scheme  as 
described  in  Definition  3 -3 ;  for 

|z|  =  1  ,  z  ^  1  ,  ita(z)  =  1  , 

(5 .11)  is  implied  by  (3.7)}  and  finally,  by  Lemma  5«1>  k  =  1  is 
never  an  inner  root  for  z  ~  1.  Thus,  (5.11)  is  verified  and  the 
theorem  is  proven. 

We  conclude  the  paper  by  proving  Lemmas  3*1  and  3*2. 

Proof  of  Lemma  3.1.  (i)  Let  the  boundary  scheme  (3.1)  be 

two-level  and  accurate  of  order  zero  at  least.  By  zero-order 
accuracy,  the  boundary  coefficients  satisfy 

m  a  m 

^  cj(-l)  “  ^  ^  * 

j=0  <j=0  j=0  ** 

hence 

(5-32)  R^z»,{^z=*«=l  =  0  * 

Since  the  boundary  scheme  is  two-level,  R(z,k  =1)  is  a  first 
degree  polynomial  in  z"^  whose  single  root,  by  (5.12),  is 
z”1  =  1.  Consequently, 

R(z,k  =  1)  {  0  ,  |  z|  =  1  ,  |z|  ^  1  , 

and  (3>7)  holds. 

(ii)  In  the  three-level  case,  R(z,k  =  1)  is  a  2nd  degree 

polynomial  in  z"1  with  real  coefficients.  By  (5;  12)  again, 

-1 

z  *  1  is  one  of  the  roots,  so  the  other  is  real  as  well, 


namely 
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(5*13)  R(z,tc  «  1)  jf  0  ,  )z|  *  1  ,  z  /  +  1  . 

Combining  (5*13)  with  our  hypothesis 

R(z  -  -1,  k  -  1)  /  0  , 

(3*7)  follows  and  the  proof  is  complete. 

Proof  of  Lemma  3.2.  As  in  the  proof  of  Lemma  5*2,  apply 
Lemma  4.1  to  the  boundary  scheme  (3-8).  Since  the  boundary  scheme 
is  right  sided  with  co(_i)  ^  0  ,  we  find  that  it  is  solvable  if 

(a)  the  difference  equations 

(5.14)  T_x  wv  S  Cj(.1)  wv+j  =  0  ,  V  =  0,1,2,..., 

have  no  nontrivial  solution  w  e  ig(x),  and 

(b)  the  hypothesis  of  the  present  lemma  is  fulfilled,  i.e., 

(5-15)  T^k)  =  cj(  ;L)  k*  +  0  ,  0  <  |k |  <  1 

Now,  it  is  well  known  that  the  most  general  solution  of  (5*14)  in 
/g(x)  is  a  combination  of  powers  of  the  roots  of  T_^(k  )  which  lie 
inside  the  unit  disc.  Thus,  if  (5*15)  holds  then  the  only  solution  of 
(5«l4)  is  the  trivial  one,  namely  (b)  implies  (a)  and  the  lemma 
follows . 
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